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Abstract 

We argue through a combination of slave boson mean field theory and the Bethe ansatz that 
the ground state of closely spaced double quantum dots in parallel coupled to a single effective 
channel are Fermi liquids. We do so by studying the dots' conductance, impurity entropy, and 
spin correlation. In particularly we find that the zero temperature conductance is characterized by 
the Friedel sum rule, a hallmark of Fermi liquid physics, and that the impurity entropy vanishes 
in the limit of zero temperature, indicating the ground state is a singlet. This conclusion is in 
opposition to a number of numerical renormalization group studies. We suggest a possible reason 
for the discrepancy. 
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I. INTRODUCTION 



Quantum dots provide a means to realize strongly correlated physics in a controlled 
setting. Because of the ability to adjust gate voltages which control both the tunnelling 
amplitudes between the dots and the connecting leads and the dots' chemical potential, 
quantum dots can be tuned to particular physical regimes. One celebrated example of said 
tuning was the first realization of Kondo physics in a single quantum dot,-i^ obtained by 
adjusting the chemical potential of the dot such that the dot was occupied by one electron. 

More generally, engineered multi-dot systems offer the ability to realize more exotic forms 
of Kondo physics. This was seen, for example, in the realization of the unstable fixed point 
of two-channel overscreened Kondo physics in a multi-dot system.- There has thus been 
considerable theoretical interest in double dots systems in different geometries, both in series 
(for example Refs. IsMlOl) and in parallel (for example Refs. Illl-l24j). In this article we focus 
on the strongly correlated physics present on the latter geometry, in particular when there 
is no direct tunneling between the dots and the dots are not capacitively coupled. Such dot 



systems have been experimentally rea^ 
of double dots in parallel see Refs. |28 



ized in numerous instances^^ — (for other realizations 



291 1. Although the dots are not directly coupled, they 



are coupled through an effective Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction. It is 
aim of this work to explore the nature of the RKKY interaction in parallel double quantum 
dots. 

A straightforward application of the RKKY interaction as typically understood would 
lead one to believe that in a closely spaced double quantum dot with two electrons present, 
one on each dot, the RKKY interaction should lead to an effective ferromagnetic coupling 
between the dots. How should this coupling reveal itself in a transport experiment, the 
typical probe of a quantum dot system? If a ferromagnetic coupling is present, one expects 
the electrons on the two dots to bind into a spin-1 impurity. If the dots are coupled to a 
single effective lead, we obtain, in effect, an underscreened Kondo effect. As the temperature 
is lowered, the single effective lead will partially screen the spin-1 impurity to an effectively 
uncoupled spin-1/2 impurity. The ground state of the system will then be a non- Fermi 
liquid doublet. In particular at small temperatures and voltages, the conductance through 
the dot will be characterized by logarithms.— This scenario has been put forth in a number 
of NRG studiesi^>i^>i^>i^ and is implicit in a number slave boson studiesi^"— of double dots 
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in parallel. 

We present contrary evidence here that this scenario is not in fact applicable at least 
for temperatures below the Kondo temperature. We argue that the ground state of closely 
spaced double dots is instead a Fermi liquid singlet. These findings are consistent with those 
of Ref. y. We do so using both the Bethe ansatz and slave boson mean field theory. It has 
been shown^iii^ that under certain conditions double dots in parallel admit an exact solution 
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using the Bethe ansatz. This exact solution, following the approach introduced in Ref. 
can be exploited to compute transport properties. In particular, the zero temperature linear 
response conductance can be computed exactly. Double quantum dots, however, only admit 
an exact solution provided their parameters satisfy certain constraints. To ensure that 
our finding of Fermi liquid behaviour is not an artifact of these constraints, we also study 
the parallel dot system using slave boson mean field theory. This allows one to study the 
sensitivity of our results to adding a second weakly coupled channel and to compute such 
quantities as the spin-spin correlation function, an object not directly computable in the 
Bethe ansatz. 

The paper is organized as follows. In Section II we detail the double dot model that 
we are interested in studying together with the approaches (Bethe ansatz and slave boson 
mean field theory) that we employ in studying this system. In Section III we present results 
on the linear response conductance through the dots both at zero temperature and finite 
temperature. We show the zero temperature conductance obeys the Friedel sum rule, a hall 
mark of Fermi liquid physics. We also study the impurity entropy at finite temperature 
showing that it vanishes in the zero temperature limit indicating the presence of a singlet. 
Finally in this section, we present results (using slave boson mean field theory alone) of 
the spin-spin correlation function. Lastly, in Section IV, we discuss the implications of our 
results and suggest a way they can be reconciled with the conflicting NRG studies. 

II. MODEL STUDIED 

We study a set of two dots arranged in parallel with two leads. The Hamiltonian for this 
system is given by 



n = -i 

lo 



/oo 
dxc\^d^ci„ + ^ Vic,{c\^^d^c + h.c) 
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FIG. 1: A schematic of the double dot system. 



+ ^^danaa + ^Uan^ania- (2.1) 
(TO a 

The Qo- specify electrons with spin a hving on the two leads, / = 1,2. The d^a specify 
electrons found on the two dots a = 1,2. Electrons can hop from the leads to dots with 
tunneling strength Via- The strength of the Coulomb repulsion on the two dots is given by 
Ua- We suppose that there is no interdot Coulomb repulsion and that tunneling between 
the two dots is negligible. A schematic of this Hamiltonian for the two dots is given in Fig. 
1. 



A. Bethe Ansatz 

The above Hamiltonian can be solved exactly via Bethe ansatz under certain conditions. 
The set of constraints that we will be particularly interested in are as follows: 

Vla/V2a = Via'/V2a'', 

Ua + 2eda = Ua' + 2eda'. (2.2) 

The first of these conditions results in only a single effective channel coupling to the two 
leads. This occurs automatically when the dot hoppings are chosen to be symmetric and 
so is commonly found in the literature.— '^'i^'^S"—'^ The second condition tells us that 
with f/i2 fixed, e^i — €^2, is also fixed. We thus must move €^1^2 in unison in order to 
maintain integrability. The final condition tells us that y/UiTi, the bare scale governing 
charge fiuctuations on the dots must be the same on both dots. 
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To solve this Hamiltonian we implement a map to even and odd channels, Ce/o = 
(yi/2,aCi ± V2/i^aC2)/ V^^a wheie = {V^^ + ^2^)72. Under the map, the Hamiltonian 
factorizes into even and odd sectors: 

/oo 
dx cl^d^Cea + V^icl^adaa + h.c) 

La aa 

aa a 

POO 

Uo = dx cl^d^Coa, (2.3) 



J o — oo 
la 

where, as can be seen, the odd sector decouples from the double dot. Using the Bethe 
ansata^*^ one can construct N-particle wave functions in the non-trivial even sector. These 
wavefunctions are characterized by N-momenta and M quantum numbers {^a\^a=\- 

The number of A^'s mark the spin quantum number of the wave function: Sz = N — 2M. 
Together the A^'s and gj's satisfy the following set of constraints: 

iq,L+i5{qj) ^ TT 9{(lj) - K+ V2 _ 

n Ag - 9{qj) + V2 ^ -TT Aa - A^ + z 
^^^Xa-9{q,)-^/2 Ha^-A.-z' 

where g{q) = ''^~'^2r~[/°^^^ • These equations are identical to the set of constraints for a single 
dot^"^ but for the form of the scattering phase S{q): 

%) = -2tan-i(5^^^^). (2.5) 

We will focus in this paper on computing transport properties in linear response. At zero 
temperature we can use the Bethe ansatz to access such transport quantities exactly.— 
We will also use the Bethe ansatz to obtain an excellent quantitatively accurate prediction 
in comparison with NRG) for the finite temperature linear response conductance (see Refs. 
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for the Bethe ansatz computation of the finite temperature linear response conductance 



46|). 



for a single dot and for the comparable NRG, Ref. 

The Bethe ansatz can be exploited to develop certain approximations that allow 
one to compute certain non-equilibrium quantities, in particular, the out-of-equilibrium 
conductance^ and the noise^ in the presence of a magnetic field. In order to obtain at 
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least qualitatively accuracy, the presence of a magnetic field is a necessity. With a magnetic 
field the Bethe ansatz for a single dot correctly predicts such features as the positioning 
of the peak in the differential magnetoconductance^ as say measured in carbon nanotube 
quantum dots.- In the absence of a magnetic field, the Bethe ansatz inspired approximation 
breaks down however.— We will, again, not consider the double dots out-of-equilibrium in 
this work. 

B. Slave Boson Mean Field Theory 

We also study the Hamiltonian fl2.ip using a slave-boson mean field theory, a well-known 
technique, applicable at sufficiently low temperatures.— The starting point is the same 
Hamiltonian f l2.ip and we will study here the Ua = oo case. The constraint of preventing 
double-occupancy on the dots is fulfilled by introducing two Lagrange multipliers Ai and A2. 
The slave boson formalism consists of writing the impurity fermionic operator on each dot 
as a combination of a pseudofermion and a boson operator: daa = bj^fcra- Here f^a is the 
pseudofermion which annihilates one "occupied state" on dot a and fej^ is a bosonic operator 
which creates an empty state on dot a. The mean field approximation consists of replacing 
the bosonic operator by its expectation value: bl^ — )■ ^6^^) = r^. Thus and Aq, together 
form four parameters which need to be determined using mean field equations. Under the 
slave boson formalism combined with the mean field approximation, Eq. (12. ip reads 



with Via = ''^aVia and ida = ^da + ^Aq. The mean field equations are the constraints for the 
dot a = 1, 2, 





(2.6) 




(2.7) 



cr 



and the equation of motion (EOM) of the boson fields. 





The above equations can be understood as arising from the conditions 

d (Hsbmft) 



6 



1 

0.8 
^0.6 
O 0-4 

0.2 




T=0 Conductance: Symmetric Dot-Lead Couplings 

Slave Bosons Bethe Ansatz 

1—^ 



1 ' 1 ' 1 ' 1 ' 




1 












1 . l/ 1 , 1 . 




-12 -10 -8 -6 


-4 -2 




-10 -8 



FIG. 2: The zero temperature conductance of a symmetrically coupled double dot computed using 
slave boson mean field theory and the Bethe ansatz. 



^ = 0. (2.9) 

Thus the reality condition on Eqn. 2.8 arises from the reality of the hopping term in the 
Hamiltonian.— For any given set of bare parameters {edaai ^a) we can compute the renor- 
malized energy {Idao) and renormalized hybridization {Via) by solving the four equations, 
Eqns. 12.71 and 12.81 While these results are mean field, they allow one to span a wide pa- 
rameter space not constrained by the requirements of integrability. For instance, we study 
asymmetrically coupled dots where two channels couple to the dot, a case not solvable by 
the Bethe Ansatz. SBMFT allows one also, unlike the Bethe ansatz, to readily study such 
quantities as the spin-spin correlation function. 

III. RESULTS 

In this section we present a number of measures as computed using both the Bethe ansatz 
and slave boson mean field theory that are indicative of the ground state of the double dot 
plus lead system. We will argue that these are consistent with the ground state of the dot 
being a singlet state, not a doublet. 
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T=0 Displaced Electrons: Symmetric Dot-Lead Couplings 
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FIG. 3: The number of electrons displaced by the dots, rid of a symmetrically coupled double dot 
computed both using SBMFT and the Bethe ansatz. In the case of SBMFT, rid is simply the dot 
occupation, ridot = ^^iai^a^i^)- ^^^^ °^ Bethe ansatz, the quantity plotted is equal to 

the dot occupation plus the 1/L correction to the electron density in the leads induced by coupling 
the dots to the leads. 



A. Zero Temperature Conductance 

The first measure we examine is the zero temperature linear response conductance, G. 
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and 



12. For the 



For the BA, G is computed as is discussed in great detail in Refs. 
SBMFT, G is computed by first solving for the variational parameters, Tq,, and ida, a = 1,2, 
and then determining the corresponding transmission amplitude via solving a one-particle 
Schrodinger equation. This procedure is detailed in Appendix B. 

If the double dot is in a singlet state, we expect G to vanish as e^i 2 are lowered, moving 
the dot into the Kondo regime. This is consistent with understanding the singlet state as a 
Fermi liquid state. If a Fermi liquid, G will obey the Friedel sum rule:— 

h 



G = 22 -7- sin^(™df7), (3.1) 



where Uda is the number of electrons of spin species a displaced by the dot. Deep in the 
Kondo regime, there will be two electrons sitting on the two dots, one of each spin species, 
i.e. rida = 1, and so G correspondingly vanishes. 
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Plotted in Fig. 2 is the zero temperature conductance as computed with BA and SBMFT 
as a function of e^i and with €di — ed2 fixed. For each computational methodology we present 
results for both e^i — ^ Fi 2 and e^i — 6^2 ^ ri,2- We see that in all cases that as e^i 2 
is lowered, the conductance vanishes. We do note however that the overall structure of the 
conductance differs as computed between the BA and SBMFT, that is to say, the SBMFT 
fails in general to describe the correct behaviour. In particular it fails to describe the 
intermediate valence regime when the distance separating the chemical potential of the two 
dots is large, i.e e^i — 6^2 ^ ri^2 • In the intermediate valence regime the conductance of the 
double dot as computed with the BA undergoes rapid changes, a consequence of interference 



between electrons tunneling off and on the dot (see Ref. Il2l ). This is not mirrored in the 
SBMFT computations which remain comparatively structureless. 

The failure of SBMFT to accurately capture the physics in the intermediate valence 
regime is seen in a related quantity, the number of electrons displaced by the dot. is the 
sum of two terms: 



cri al 



(3.2) 



The first term is simply the occupancy of the double dots while the second term measures the 
deviation of electron density in the leads due to coupling the dots and the leads. In SBMFT 
this term is zero due to the mean field nature of its approximation. However in BA this term 
is non-zero. While we cannot compute it directly, the BA gives us the ability to compute 
Ud- And as plotted on the r.h.s. of Fig. 3, we see that can be negative. As ^^{d^^dai) 
is manifestly a positive quantity, we know that as computed by the BA, the second term in 
Eq.f i3.2p is non-zero and in fact is negative (at least in the case e^i — ed2 ^ ri,2)- From Fig. 
3 we see however that Ud for both SBMFT and BA tends to two as the system enters the 
Kondo regime (where two electrons sit on the two dots). 

One advantage the SBMFT does offer over the BA is that it allows us to compute trans- 
port quantities beyond the integrable parameter regime delineated in Eq. (12. 2p . It was 
argued in RefJl2l that small deviations away from this parameter space should not qualita- 
tively change transport properties. In Fig. 4 we test this in the Kondo regime (where we 
expect SBMFT to be at its most accurate) computing the conductance while adjusting the 
dot-lead hopping parameters in such a way that we move from a case where only one effective 
channel couples to the quantum dot (i.e. Vij = 1) to a case where two channels couple to 
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T=0 Conductance of Non- Symmetrically Coupled Dots 
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FIG. 4: The zero temperature conductance of asymmetrically coupled double dot computed using 
slave boson mean field theory. The conductance is plotted as a function of Vu. The remaining 
dot-lead hopping strengths are all set to 1 while e^i = —4.7 and €^2 = —4.6. 

the dot (Vii > 1 and V12 = V21 = V22 = !)• We see when the second channel is only weakly 
coupled to the dot, the conductance remains near its one-channel value, i.e. G = Oe'^/h. 
Only once Vu appreciably deviates from 1 do we see a corresponding deviation in G. We 
note that this continuous behaviour is also consistent with the one-channel dot-lead ground 
state being a singlet. If it were instead a doublet, coupling a second channel into the system 
would lead to a discontinuous change as the second channel, no matter how weakly coupled, 
would immediately screen the free effective spin-1/2. 

Finally in this section we consider the behaviour of the conductance and displaced elec- 
trons when edi = ed2- We see from Fig. 5 that the same qualitative behaviour in both 
quantities is found using the Bethe ansatz and using SBMFT. Namely, the displaced elec- 
tron number tends to 2 while the conductance G tends to zero as e^i = €^2 goes to zero. 
While these measures are the same in the two computational methods, there is a quantity 
that sharply distinguishes the two (and so shows a failure of SBMFT even in the Kondo 
regime at zero temperature). This quantity is the low lying density of states on the dots, 
Pdifjj)- For the case of e^i = 6^2, the BA shows that pd{<^) for u on the order of the Kondo 
temperature, vanishes.— However the SBMFT shows that at this energy scale there 
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FIG. 5: The conductance and the number of displaced electrons as a function of edi{=ed2) as 
computed using SBMFT and the Bethe ansatz. 



exists non-negligible spectral weight. In Fig. 6, we plot Pd{uj) as defined by 

retarded' 

icr 

The agreement on the qualitative features of Ud and G between the two methodologies is 
then a coincidence (to a degree). In both cases the ground state is Fermi liquid like and so 
G follows the Friedel sum rule which necessarily mandates that the conductance vanish with 
two electrons on the two dots. But the nature of the Fermi liquid in each case as predicted by 
the methodologies is much different. SBMFT predicts the scale of the low lying excitations 
is Tk while the BA finds that for the special case of e^i = ed2 (and only for this caso^) that 
the fundamental energy scale corresponds to the bare energies scales in the problem, i.e. U 
and F. 



B. Finite Temperature Conductance 

We now consider the finite temperature conductance. Plotted in Fig. 7 are traces for 
G{T) for double dots with both le^i — 6^2! ^ ^1,2 and \edi — 6^2! ^1,2 as computed with 
both the Bethe ansatz as well as SBMFT. For a Fermi liquid ground state we expect that 
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-2 -10 1 

FIG. 6: A plot of the low energy density of states, p{uj) for edi = €d2 in the Kondo regime as 
computed using SBMFT. As we are argue in the text, this is an artifact of SBMFT (the BA shows 
that in this case Pd{^) vanishesi^). 



at low temperatures the conductance deviate from its zero temperature value by and 
we see that behaviour in both cases. From Fig. 7 we see that in both treatments, the 
finite temperature conductance for the double dots initial rises with increasing temperature 
to an appreciable fraction of the unitary maximum and thereafter decreases in an uniform 
manner (regardless of the bare level separation). (For SBMFT we have defined the Kondo 
temperature as where this peak in G{T) occurs while for the BA the Kondo temperature we 



12.) We however 



employ the analytic expression for Tk of the double dots derived in Ref. 
also see that there are qualitative differences between SBMFT and the Bethe ansatz. The 
peak in the conductance computed using SBMFT peaks at a value far closer to the unitary 
maximum than does the Bethe ansatz. And we also see that the conductance as computed 
in the SBMFT drops off far more rapidly than it does in the BA (particularly at large level 
separation). We however believe this is unphysical and akin to the pathologies that SBMFT 
is known to exhibit at higher temperatures and energy scales.—"— 



As was demonstrated in Ref. 
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the conductance at finite but small T is quadratic in 
T while at large T the conductance is logarithmic (going as 1/ log^(T/Tx)). The peak in 
conductance at finite T is then a result of the conductance vanishing in the low and high 
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temperature limits. These conductance profiles are similar to the those predicted in Ref. 
for multi-dots coupled to two electron channels. However the physics there is much different: 
the non-monotonicity in G{T) predicted in Ref. |4l| is due to the presence of the two channels 
and because they couple to the dots with different strengths, they screen a 5* > 1/2 state in 
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FIG. 7: The linear response conductance as a function of temperature of a symmetrically coupled 
double dot computed using both slave boson mean field theory and the Bethe ansatz. 



stages. 



C. Spin-Spin Correlation Function 

We present the static spin-spin correlation function as a function of e^i and as computed 
using SBMFT in Fig. 8. With two electrons on the dot the value of (5*1 ■5*2) can vary between 
—3/4 (if these two electrons are bound in a singlet state) to 1/4 (if the two electrons find 
themselves in a triplet state). We see in Fig. 8 that generically the value of the correlation 
function in the Kondo regime (for the relevant values of e^i see Fig. 2) that (5*1 • S2) tends 
to 0. This however should not necessarily be interpreted as the dots being closer to a triplet 
state than a singlet state. In determining the overall state of the system, (5*1 ■ S2), is not 
necessarily a good measure. We can see this by considering a simple toy example. 

Imagine a system of four spins, two associated with the dots, | t)di and | t)d2, and two 
associated with leads | and | t)i2- And first suppose the system is in a singlet state. 
Two ways this singlet state can be formed are 

Isinglet 1) = i(| t)<ii| i)d2 - I t)d2) ® (I t)zi| i)i2 - I t);2); 
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FIG. 8: The spin-spin correlation function of symmetrically coupled double dot computed using 
slave boson mean field theory. 



Isinglet 2) = -(I - I t)di) ® (I t)i2\ i)d2 - I i)i2\ t)d2). (3.3) 

We see that the expectation value, (5*1 ■ S2), of these two states is considerably different: 

(singlet ll^i ■ S'2 1 singlet 1) = -3/4 



(singlet 2151 • 5*2 1 singlet 2) = 0. 



(3.4) 



Now suppose the system is in a triplet state and suppose its projection is 1. Again 
there are two inequivalent ways this state can be formed: 

Itriplet 1) = -^1 t)di\ t)d2 ® (I t)n| 1)^2 - | i)zi| t);2); 



Itriplet 2) = —I t)n| t)di ® (| t)^2| i)d2 - \ i)i2\ t)d2). 



The expectation of these two values is 



(triplet l|5i • ^aj triplet 1) = 1/4 



(3.5) 



(triplet 2\Si ■ ^altriplet 2) = 0. 



(3.6) 
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FIG. 9: The impurity entropy as a function of temperature of a symmetrically coupled double dot 
computed using slave boson mean field theory and the Bethe ansatz. 



We thus see that when the system's state is such that (5*1 ■ = 0, it can be either a singlet 
or a triplet equally. We thus end with a more reliable measure of the dot's internal degrees 
of freedom: the impurity entropy. 



D. Impurity Entropy 

The final set of computations we present in this section are of the impurity entropy of 
the double dots. In Fig. 9 we plot results coming from SBMFT and the BA for both large 
and small level separation. (The derivation of the impurity entropy in the context of the BA 
is found in Appendix A.) We see that in all cases the impurity entropy vanishes as T — > 0. 
This then implies the ground state of the double dot system is a singlet. If it were a triplet 
state, the T — )■ limit would lead to Simp = log(2). 



IV. DISCUSSION AND CONCLUSIONS 



We have presented a number of arguments that the ground state of a double dot near 
the particle-hole symmetric point (i.e. when there are nearly two electrons on the dots) is 
in a singlet Fermi liquid state. In particular we have shown that the conductance in this 
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limit vanishes, in accordance with the Friedel sum rule and that the impurity entropy also 
vanishes, in agreement with the ground state being a singlet. 

These conclusions however disagree with a number of NRG studies. In Refs. 



13 



and 



14 



it is found that a double quantum dot (with the dots closely spaced) carrying two electrons 
is in a spin-triplet state, has a conductance corresponding to the unitary maximurn, anc 
is correspondingly a non- Fermi liquid. Similar conclusions are reached in Refs. [igIIi? . 
The basic rational invoked for observing this physics is that with dots closely spaced, a 
ferromagnetic RKKY interaction is present which binds the two electrons on the dot into a 
spin triplet. Consequently the ground state of the double dot is that of an underscreened 
spin-1 Kondo impurity. 

We offer a possible reason for the discrepancy that we find with these studies. The 
Anderson Hamiltonian typically considered in these studies is of the form: 

H = Hlead + Hdot + Hiead-dot'-i 



Hlead ^ ^ J kdkal^ak^; 



i=l,2\a 



Hdot-iead = D'^^J2 f dk V^idlttka + (4.1) 

where we are using the conventions of Ref . |48| in writing down the NRG Hamiltonian. Before 
implementing the NRG algorithm, one adopts a logarithmic basis for the lead electrons, 

aka = ^anpa'iPnpik) + bnpa'iPnpik), (4-2) 
np 

where 

^i(^) = ^^'^ <±^<A , ^^_3^ 



if k is not within the above interval, 



and Wn is given by 

Here A is a parameter less than one. This change of basis transforms Hiead and Hiead-dot 
into 

Hlead = y (1 + A"^) A""(«npa«npfT " b^p^hnpa) 

np 
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D{l-A-^)^^ .27n{p' -p).. . . ^ . 

2m ^ ^ exp(^-^^)(aTp,a„pv - bl^^bnp'a)-—; 

Hi,aci-,ot = D^'\l - A-i)i/2 5^ A-"/V.((4o. + + h.c). (4.5) 

incF 

Typically the approximation that is now made in most NRG treatments (and seems to 
have been made in the above references) is that the p 7^ modes of the logarithmic basis 
are dropped, not least because these modes do not directly couple to dot. In the single 
dot case, where this approximation was first made,— 1^ this was found to be a reasonable 
approximation. However in the double dot case it is not a priori obvious that this is the 
case. In particular if one finds an underscreened Kondo effect one might ask whether that 
additional modes {p 7^ 0) might serve to provide additional screening channels. 

Let us then consider the NRG Hamiltonian that would arise if both the p = and p = ±1 
modes were kept, if/ead can then be trivially diagonalized using the combination 

= —j= [5ania + (2 =F 6i)e^a„o<7 - (4 ± 3i)e^^an-ia] , 
3v 10 

with 9 = — ^^^-1 . Note that corresponding transformation {s ^ 6} is omitted for brevity. 
This yields 



H -DVA-4l±^rtr , 27r(l + A) + 3(1 - A) 



27r(l + A) -3(1-A)^^ 



47rA 



Hdot-lead — Vr 



m 

ina 



AttA 

The corresponding transformation of dot-lead Hamiltonian is 

\/f ~ ^) ^ ~ \^Oadia + {r ^ s} + h.c 

(4.6) 

with Vni = D^1- A-i)U-5e^yi. 

We see upon this diagonalization that three channels of electrons, t„{o,±},<t and Sn{o,±}a, 
couple to the dot. And because of the nature of the logarithmic basis, we see that +, — 
and variables can be arbitrarily close to the Fermi surface (due to the presence of the A"" 
factor) and so should all contribute to Kondo screening. 
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This analysis suggests that in a situation with two electrons on the double dots which 
are bound into a triplet by a putative RKKY coupling, there are nonetheless (at least) three 
available screening channels, at least in the NRG reduction of the Anderson Hamiltonian. 
And so the problem would seem to be not one of an underscreened Kondo effect but an 
exactly screened Kondo effect. We do however note that while viewing the double dot 
system with two dot electrons as an exactly screened Kondo effect is consistent with our 
Fermi liquid findings, it is not clear under what conditions it would be correct to think of 
the two electron on the dots as ever forming a triplet. While any perturbatively generated 
RKKY interaction will generically be larger than the exponentially small Kondo screening 
scale, there is a question of whether the zero temperature perturbation theory underlying 
any RKKY estimate is convergent because the system is gapless.— 

While we have focused here on NRG treatments of double dots close to their particle-hole 
symmetric point, similar analyses of any multi-dot system might suggest that whenever the 
dot degrees of freedom exceed S = 1/2, it is not possible to ignore the p ^ modes of the 
logarithmic basis. 



Acknowledgements: MK was supported by the NSF under grant no. DMR-0906866. 
RMK acknowledges support by the US DOE under contract number DE-AC02-98 CH 10886. 
We thank Chung-Hou Chung, Alexei Tsvelik, Piers Coleman and Leonid Glazman for useful 
discussions. 



Appendix A: Analysis of the Ground State Entropy via the TBA Equations 



Here it is demonstrated that the ground state entropy of the double dot system at zero 
temperature is zero and thus the ground state is a singlet. The procedure outlined below 



42 



can be found for a single dot in Section 8.3.3 of Ref. 

We start with the observation that the free energy of the system can be expressed as 
sums over all excitations in the system, that is, over all possible solutions of the Bethe 



ansatz equations (see Eqn. 10 of Ref. Il2[ ). Specifically it takes the form 



E-TS, 



(Al) 
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where the energy of the system equals 

oo „ 

E = I dkp{k)k + J2 / rfAa;(A)eo„(A), 



n=l 



(A2) 



where 



eon(A) = -n{2edi + Ui/2) + 2 / a^iX - g{k))g{ky, 
an{x) 



9{k) 



Tim? + ' 
{k - e,i - ^tf 



(A3) 



and the entropy, S", is given by 



S 



dk 



(Ak) + m) log(p(»:) + pm - p(k) \ogp(k) - ^(A^) \og(p(k)) 



+E 

n=0 

oo 



n=0 



(cr„(A) + (T„(A)) log(cr„(A) + cr„(A)) - a-„(A) logcr„(A) - o-,i(A) log((T,i(A)) 
«(A) + a;(A)) log(a;(A) + d'M) - ^n(A) loga^A) - ^^(A) log(a;(A)) 



■ (A4) 



Here p(/c), (T„(A), and o"^(A) are the particle densities while 5"„(A), and cr^(A) are the 
hole densities of the various excitations (i.e. solutions of the Bethe ansatz equations). The 
particle and hole densities can be shown to obey the following equations: 



p{k)+P{k) = Po{k)-g'{k) dXs{X-g{k)){a,{X) + a[{X)), 



(A5) 



where 



Poik) 
A{k) 
A„(A) 



^ + + g'ik) I dXsiX - gik))i-^e',,iX) + ^Ai(A)); 



2tt 



dkS{k); 



1 



--ReA{V2U^,{X + ^)y, 



TT 



Six] 



1 

2 cosh(7rx) 



(A6) 
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and 



dnW + MX) = j rfA's(A-A'))(a„+i(A') + an-i(A')) + 5„i j dkp{k)s{X - g{k))- 

^nW + a'SX) = I dX's{X - X'))K^,{X') + ~a'^_,{X')) + D„(A), (A7) 
where 

Dn{X) = 5„i ! dkp{k)s{X - g{k)) - [ dX's{X - A'))(A„+i(y) + A„„i(A')); (A8) 



One sees that the density equations have source terms that involve a bulk piece and a piece 
scaling as 1/L, where L is the system size. 

One defines the energies of the excitations at finite temperature via the relations, 

<k)=T\og{^-^); e„(A)=Tlog(^); <(A) = T log(^). (A9) 

These energies are given by the relations e(/c),en(A), and e^(A) which are governed by the 
equations, 

eik) = k + I dXeoiiX)s{X-gik)) + T J rfAs(A - ^(A:)) log ; 
e„(A) = 5„ir J dkg'ik)siX-gik))\ogini-e{k))) 

-T I dA'.(A-A')log(n(6Li(A))n(6'„+,(A))); 
e:,(A) = 6n,T I dkg\k)s{X- g{k))\og{n{e{k))) 

-T [ dX's{X-X')log{n{e'^_,{X))n{e'{X))), (AlO) 



where n{x) = (1 + exp(x/T))~^ is the Fermi function. The equations for the e's and the 
bulk pieces of the densities (i.e. the pieces not scaling as 1/L) are the same as for a single 
level dot. As noted in the manuscript, the Bethe ansatz equations for the double dots in 
parallel are identical to the single level dot up to the impurity scattering phase. 

Substituting the energies and densities in the expression for the free energy, one can 
rewrite it in a much more simple fashion: 

dkpo{k)\og{n{-e{k)))+T / rfA p^,(A) log(r2(6'„(A))); 

^ n=0 
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/Hk 1 
-s{X-gik)) + -DM- (All) 

Here Egg is the ground state energy of the system. 

One now wants to consider how Q behaves at T — 0. If one can show that the leading 
order correction in Q at low temperatures is then as S = —dr^ one will have shown that 
the entropy vanishes as T — )■ 0, and so the ground state of the system is a singlet. 

In order to see that Q has no term linear in T, it is sufficient to consider the zero 
temperature values of the energies, e{k) and e^(A). At the particle-hole symmetric point, 
one has 

e{k,T = 0) > , for all k; 



e[{X,T = 0) < , for all A; 



e;(A,T = 0) = , n > 1, for all A. (A12) 

If one substitutes these expressions into the expression for Q and uses the fact that 
/ dXpl,{X) = 0,n>l, one sees that Q = Egs + 0{T^). 

Now if one is away from the particle-hole symmetric point, one has instead 

e{k,T = 0) > , for all k; 

e;(A, T = 0) = n(y + ei) , for all A. (A13) 
Now while e'i(A) is neither solely positive nor solely negative at zero temperature, its leading 



order finite temperature correction is (see Section 8.3.7 of Ref. |42| ) 



e[{X,T) = e[{X,T = 0) + O{T^). (A14) 

Substituting these forms of the energies into the expression for the free energy, one again 
sees that there is no term in Q that is linear in T. 
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Appendix B: Derivation of the Conductance in SBMFT 



Here we present a derivation of the conductance in the general case of asymmetrically 
coupled dots. To determine the conductance we solve the one-particle Schrodinger equation 
of the SBMFT Hamiltonian, Hsbmft\''P) = E\iIj) where lijj) equals 



+00 



> = I dxgi{x)c\{x)\Q > + / dxg2{x)c\{x)\Q > 

+ei4|0 > +e24|0 > • 

This gives the following four equations: 

- id^gi{x) + eiVii5{x) + 62^12^(3;) = Egi{x); 

-id.j,g2{x) + 61^^21^(3;) + 62^^225(3;) = Eg2{x); 

{id, -E)e, + ViigiiO) + ^1^72(0) = 0; 

{id, -E)e2 + ^22^?2(0) + V^i2^7i(0) = 0. 

We then take the functions 5'i,2(a;) found in the one particle wavefunction 
following form: 

g,{x) = e'^''{e{-x)+R^^9{x)); 
e^^^^(a;)Ti2. 



92[X) 



(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 

to be of the 

(B6) 
(B7) 



Substituting the ansatz (Eqns. IB6I and Eq. IB7p into the above four equations, we obtain 
four equations from which one can solve for T12. The conductance G, equal to G = ^Ti2Tj*2, 



is then 



G 



2^D 
h N 



N = 16 



^12^22edi + ViiV2lid2 



D = 8 



V'llV^12 + V^2lV^22 



id,id2 + 16ed,4 



V12V21 — V11V22 



+4e 



~2 
di 



v,\ + V2I 
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y2 + T/2 
^11 ^ ^21 



We have also computed the transmission amplitude using Ref. 



T = Tr <^ CTrG'Tl 



50 



(B8) 



(B9) 
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where G"/^ are advanced and retarded Green's function matrix and Tr and Tl are defined 
by 

V21 V22 



R 



Vo,V, 



21 ^^22 



'^22 



Vi,V,2 Vi 



12 



(BIO) 



We find that this trace formula (Eq. IB9I) gives exactly the same result. 

In the symmetric case (i.e. e^i = ed2 = Cd), the above mean field equations reduce to two 
equations given by: 

^2f 



r 1 1 
= - - 1 + -tan"^ 

r 71 



id - ed 1 
= — - — + - log 

i 71 



(Bll) 



Here one can compute by solving Tk = ed — [^] which in the limit <^ \ed\ 

gives Tfc = DeTvr. In this case the Kondo limit is reached when the gate voltage reaches 
ed = -4.45 (with T = 1 and D = 100). 
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